From superoperator formalism to nonequilibrium Thermo 

Field Dynamics 



Y. Nakamura a , Y. Yamanaka a 

" Department of Electronic and Photonic Systems, Waseda University, Tokyo 169-8555, Japan 



Abstract 

Emphasizing that the specification of the representation space or the quasiparticle picture 
is essential in nonequilibrium quantum field system, we have constructed the unique 
unperturbed representation of the interaction picture in the superoperator formalism. To 
achieve it, we put the three basic requirements (the existence of the quasiparticle picture 
at each instant of time, the macroscopic causality and the relaxation to equilibrium). 
From the resultant representation follows the formulation of nonequilibrium Thermo 
' Field Dynamics (TFD). The two parameters, the number distribution and excitation 

energy, characterizing the representation, are to be determined by the renormalization 
condition. While we point out that the diagonalization condition by Chu and Umezawa 
is inconsistent with the equilibrium theory, we propose a new renormalization condition 
as a generalization of the on-shell renormalization on the self-energy which derives the 
quantum transport equation and determines the renormalized excitation energy. 
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1. Introduction 

A choice of the representation space, or the Fock space, is essential in quantum field 
theory. Without specifying it , it would be impossible to calculate matrix elements of 
field operators, namely physical quantities. In quantum field theory of vacuum, we can 
rely on the picture of the stable asymptotic particle on the stable and unique vacuum, 
and the representation space is the Fock space associated with the asymptotic particles. 
Because the concept of the asymptotic field is achieved by the disappearance of overlap 
of the wave-packets describing the particles in the infinite past and future, it is clearly 
invalid in the thermal system. If we insist on the existence of the asymptotic particles 
in thermal system, the S- matrix is unity (so-called "trivial S-matrix" ) [l|, 2| ■ To this one 
may add the physical argument that a Hamiltonian for thermal system is not bounded 
below (see e.g. Eq. and that a stable particle picture cannot be admitted, which is 

seen from non-vanishing imaginary part of the self-energy under the thermal situation. 
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Hence every particle becomes thermally unstable and should be called a quasiparticle. 
What we aim in this paper is to construct the quasiparticle picture explicitly in the 
unperturbed representation of the interaction picture under thermal situation, which 
supplants the asymptotic particle picture in the vacuum theory. 

There are known two nonequilibrium thermal field theories, i.e., the closed time path 
(CTP) formalism Q and Thermal Field Dynamics (TFD) 0. While the CTP formal- 
ism is widely used, we have been employing the TFD formalism because the concept of 
the representation space is clear even in nonequilibrium situation. In TFD, which is a 
real-time canonical formalism of quantum field theory, the thermal degree of freedom is 
introduced through doubling each degree of freedom, and the mixed state expectation in 
the density matrix formalism is replaced with an average of a pure state vacuum, called 
the thermal vacuum. A well-defined quasiparticle picture is constructed on the doubled 
Fock space, and quasiparticle operators which diagonalize the unperturbed TFD Hamil- 
tonian is defined in a self-consistent manner. A time-dependent number distribution is 
introduced as an unknown parameter, and the self-consistent renormalization condition 
Q derives an equation for it, i.e. the quantum transport equation, which reduces to the 
well-known quantum Boltzmann equation in the Markovian limit. 

So far we have investigated the cold atomic gas systems without and with Bose- 
Einstein condensate [y, M, and have derived the quantum transport equation in 
nonequilibrium TFD [9|, [l0(. The appropriate choice of the quasiparticle picture was 
essential there. As an important result of our previous works, our transport equation 
in the presence of Bose-Einstein condensate obtains an additional collision term, called 
the triple production term, which is absent in the other approaches llj, [12 . El. While 



the triple production term vanishes for the stable condensate, it remains non-vanishing 
in case of the thermally unstable condensate to prevent the system from equilibrating. 
Such instability, called the Landau instability and caused by collisions in which the 
neg ative-energy quasiparticles participate, has been observed in the cold atomic system 
[14j. An approach using an inappropriate particle picture, for example the particle pic- 
ture of original atom with energy spectrum E(x,p) = p 2 /2m + V(x), takes account of 
no negative-energy particle and is inadequate to describe the Landau instability. This is 
an example showing the importance of the choice of the quasiparticle picture. 

In this paper we take a position that a representation space corresponding to a par- 
ticular quasiparticle picture is chosen by taking an unperturbed representation in the 
interaction picture. This is done in the superoperator formalism [15j, following from 
the density matrix formalism. Then we put the three basic requirements which thermal 
field theory with the coexistence of microscopic and time-dependent macroscopic quan- 
tities is desired to fulfill, and acquire the unique unperturbed representation in which 
the unperturbed time-dependent density matrix has the geometrical distribution as the 
equilibrium one. The discussions of this paper is confined to cases of no spontaneous 
breakdown of symmetry, i.e. without condensates. As will be easily seen, the formula- 
tion of nonequilibrium TFD is derived from this representation. It should be emphasized 
that what we have done in this paper is not a mere re-derivation of nonequilibrium TFD, 
but provides a new perspective and understanding of it. First, it is shown that the for- 
mulation of nonequilibrium TFD using the thermal vacuum and time-dependent thermal 
Bogoliubov transformation, which has been assumed without a sound rationale, corre- 
sponds to the unique consistent choice of the unperturbed representation. We remark 
that the Feynman diagram method, a very powerful tool of quantum field theory though 
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it is never guaranteed in thermal field theory for nonequilibrium system, is available ow- 
ing to it. Second, we point out an undesirable property of the renormalization condition 
proposed by Chu and Umezawa Q, called the diagonalization condition. Explicitly the 
diagonalization condition is consistent only in the leading order, since it conflicts with 
the equilibrium theory in higher orders. In the context of the discussions in the present 
paper, we find and propose a new self-consistent renormalization condition which over- 
comes the difficulty. Our renormalization condition reduces to the diagonalization one in 
the leading order, nevertheless it is consistent with the equilibrium theory even in higher 
orders. 

This paper is organized as follows. We briefly review the formulation of the super- 
operator formalism in Section 2. In Section 3, the Liouville-von Neumann equation is 
introduced in terms of superoperator formalism, and the Schrodinger, Heisenberg and 
interaction pictures are introduced in the Liouville space. In Section 4, the unperturbed 
quasiparticle representation is constructed from the three basic requirements. The unique 
structure of the unperturbed density matrix is derived from the requirements in Section 
5, which has been just an assumption in nonequilibrium TFD so far. The equivalence of 
the superoperator formalism to nonequilibrium TFD is shown in Section 6. In Section 
7, a new self-consistent renormalization condition is proposed. Section 8 is devoted to 
summary and discussions. 

2. Superoperator formalism 

We briefly review the superoperator formalism for our discussions below, following 
Schmutz 

Let a quantum system be described by a set of annihilation- and creation-operators, 
[aj , 4]a = Sjk , [aj , a k ] a = [aj , a\] a = (1) 

with 

[A, B] a = AB - aBA, 

I 1 for bosons 
a = < 

— 1 for fermions . 



(2) 



The Hilbert space (Fock space), denoted by H, is spanned by the complete orthonormal 
set of basic vectors which are generated from operations of aj on the vacuum |0) , 

|m) = |mi,ma, • ■ • ) = J] (a})™ 3 |0) , (3) 

j V 3 ■ 

(m\m'} = 8 mm > , }Jm)(m\ = I . (4) 

m 

The linear operators acting on % form a linear space, called the Liouville space % . 
We denote its element by a double- ket | A~j , corresponding to the operator A acting on 
Ti, where the notation of the double-ket |-} is used to distinguish it from the vector |-) 
of %. The inner product of two elements is defined by 

(A\B}=Tr[A*B]. (5) 
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The complete orthonormal basis of H is {|m, nj> = | \m) (n|))} , 

{m,n\m',n'} = 5 mm >5 nn > , J~] \m,n)(m,n\ = I, (6) 

where (m, n| = ((|m)(?i|| . For each operator acting A on H, one has a superstate |Aj, 



in./ 1 



(m, n | A) = (m|A|n) . 
We introduce a concept of superoperator, operating on the superstates, and a special 
set of superoperators, dj and dj, by 

a, |m, n)) = |a :) '|m)(n|} , aj|m,7j)) = \foo ii ~ v \ \m) (n|a])) , (8) 

where /x = ■ rrij and ^ = ■ n -i . T his phase convention of the fermionic tilde operator 
is different from Schmutz's one [TBI], but is adopted from the references [l6|, 0, E3) EH- 
From the Hermitian conjugates of Eq. (|5|). 

(m,n|at = (a i |m)(n|| > «w,n|at = Vo : cr^- , ' +1 (|m)(n|ot| , (9) 

and the definition of the inner-product ([3]), one consistently obtains 

fit|m,n) = |ot|m)<n|>, a]\m,n) = sftj(T>*- v \\m)(n\a i ) . (10) 

The superoperators dj and dj simulate the operations of dj on |m)(n| from the left and 
right, respectively. 

It follows from Eqs. © and ([TO]) that the superoperators obey the (anti)-commutation 
relations, 

[dj , 4]cr = [dj , 4]<T = (Sjfc , 

[6j , a fc ] CT = [a] , d\] a = [dj , d k ] a = [a] , aj.]<r = . (11) 

The basic vectors of T-L are generated from operations of the creation superoperators on 
the supervacuum 0, 0^, 

\m, n> = (^r 2 (n ^=1 (n ^) 4 • • • 4A • • • 4> o> , ^ 

and are the eigenstates of supernumbers, 

d'-dj |to, = rrij\m, ?i]> , aja^m, = 7ij|m, n]> . (13) 
The identity operator / in ^ gives the superstate |/)), 

|7»=^|m,m». (14) 



This, together with Eqs. (|8J and (fT0|) . yields the relations, 

fa - V^afj \I) = fa - VSii) I J» = . (15) 

Suppose an operator A — A(a\a) in W, then we define two superoperators by 

A = A{a\a), A = A*(a\a), (16) 

and refer to them as non-tilde and tilde operators, respectively. The relation between 
them is summarized as the tilde conjugation rules, 

(iii 2 ) " = ilia , 

(ciii + c 2 i 2 ) ~ = cjii + c* 2 A 2 , 

(it)~ = it, 

(i)~ = i, (17) 



where q are c-numbers. Note that the last property and Eq. (|15p are due to the phase 
convention for fermionic operator in Eq. (|5J| and differ from the conventional ones [J] . It 
can be shown from Eq. (fT5|) that 

m)=l A \ I X forbosonicA as) 

Iij4 7^ for fermionic A . 
The tilde conjugation may be applied to superstates, and we find 

\m,nf = {^) {p+ " ]2 \n,m) . (19) 
Therefore |/)) becomes tilde-invariant, 

|/»~=|/». (20) 

When A is bosonic, we can show 

W= (i|/»)~=i|/»=it|/» = |At>, (21) 
from Eqs. (gU]) and (fTS]). This implies that 

A* = A is equivalent to \A} — \A} for bosonic A . (22) 

3. Liouville von Neumann equation in the Schrodinger, Heisenberg and in- 
teraction pictures 

Consider a quantum system under thermal situation in the Schrodinger picture, de- 
scribed by the density matrix ps(t)- Its dynamical equation is the Liouville- von Neumann 
equation, 

ij t p s (t) = [H s ,p s (t)], (23) 
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with the original total Hamiltonian H$ which is Hermitian. We take the unit of h = 1 
throughout this paper. We convert (1231) to the equation in the superoperator formalism, 
described in the previous section: 



i-\ps(t)} = H s \p s (t)}. 



dt 

where the total super-Hamiltonian is 

Hs = Hs — H$ ■ 



(24) 



(25) 



The Heisenberg picture is constructed as usual. We introduce the time-evolution 
operator by 

(26) 



U (t, t ) = exp (-iH s (t - t ] 



and to is the coincidence-time of the various pictures. Then the time-independent state 
\ph) is related to |ps(*)} as 



\ps{t)) = U{t,to)\ PH ), 
while the Heisenberg superoperator is given by 

A H {t) = U- 1 (t,t )A s U(t,t ), 
and satisfies the Heisenberg equation, 

d 



i—A H (t) = A H (t), H(a H (t), 
at l 



(27) 



(28) 



(29) 



In order to move to the interaction picture, the unperturbed super-Hamiltonian H u 
needs to be specified. Here are the two remarks regarding to it. Firstly it generally 
differs from the free super-Hamiltonian Hq = Hq — Hq , constructed from quantum field 
theory of vacuum. Secondly it involves some time-dependent parameters, say ry(t), which 
will turn out to be the mean number rij(t) and the renormalized excitation energy uij(t), 
and therefore depends explicitly on time, in contrast to the usual unperturbed Hamilto- 
nian. This is characteristic of our approach to describe time-dependent nonequilibrium 
processes, which will be clarified in later sections. Once H u is specified, we introduce the 
time-dependent superoperator A(t) in the interaction picture by 



A{t) = UQ X {t,to)A s Uo{t,to) 



where 



i-T;U (t,to) = H u (a s , > 



r)(t))U (t,t ) , 



(30) 



(31) 



and ri(t) inside H u represents time-dependent parameters, mentioned above. The time- 
evolution operator Uo(t,to) is 



Uo(t,t ) = T 



exp 



i / H u (a s , • • • : ds 



(32) 



where the symbol T stands for a time-ordered product, 

T [A{h)B(t 2 )} = 6{h - t 2 )A(h)B(t 2 ) + a6{t 2 - ti)B(t 2 )4(ti) . (33) 
The equations for {&j(t) , aj(t) , dj(t) , a](i)} are given by 

d 



dt 



(34) 



and so on. Here the unperturbed super-Hamiltonian in the interaction picture is denoted 
simply by H u (t) , 



H u (t) = H u (a(t),--- :r)(t)) = Uo\t,t )H u (a s ,- ■ ■ : f/(*))t>b(t, *o) - 



(35) 



Note that the unperturbed Hamiltonians in both the interaction and Schrddingcr pictures 
do not coincide with each other in general due to the time-dependent parameters T](t), 



others = . 



The equal-time (anti)-commutation relations are 

4 (*) = 5j(*)>4(*) = S 3k, 

L J a L J (7 

The state in the interaction picture is 

\p(t)} = Uo\t,t )\p s (t)} = V(t,t )\ps(t )} , 
V(t,to) = Uo 1 (t,to)U(t,t ), 

ijV{t,t ) = Hi(t)V(t,to), 

with the interaction super-Hamiltonian, 

H I (t) = H(a(t),---)-H u (a(t),---:v(t))- 

Formally V(t,t ) is 



(36) 



(37) 



V(t,t ) = T 



exp 



Hi{s)da 



The causal two-point Green's function is defined in the Heiscnbcrg picture, 



G(h,t 2 ) 



-i (IT 



A H (ti)B H (t 2 ) 



\Ph] 



(38) 

(39) 
(40) 

(41) 



which is — iTr [T [Ah {t\)BH{t 2 )] pn] in the density matrix formalism, when An(t) = 
An(t) and Bn(t) = Bn(t) ■ We rewrite it in terms of the operators in the interaction 
picture, putting to = — oo, 



G(t 1 ,t 2 ) = -i(/|S- 1 T SA(ti)B(t 2 ) 
where the S-matrix superoperator is 

S = V(oo, -co) . 
7 



(42) 



(43) 



4. Unperturbed representation 

In this section we attempt to construct the unperturbed representation in the super- 
operator formalism from general requirements for quantum field system under thermal 
situations. 

We maintain that physical quantities for a quantum field system, following from 
matrix elements of the field operators, can be calculated only when the Fock space 
representing the abstract field operators is specified. In quantum field theory of vacuum, 
the existence of the stable and unique vacuum and the stable asymptotic particles makes 
the choice of the Fock space unambiguous, namely the Fock space of asymptotic fields 
is chosen. With to = 0, the vacuum in the Heisenberg picture is related to that of the 
asymptotic fields by the Gell-Mann-Low relation. Alternatively, one may take to = — oo, 
then both of the vacua coincide with each other. In addition, because of (0|>S _1 oc (0| 
(S: the S-matrix operator, |0): the asymptotic vacuum) we have the Feynman diagram 
method in calculating the causal Green's functions. 

The concepts of the vacuum as the ground state and the asymptotic particles are 
invalid in thermal situation. As is seen from Eq. the total Hamiltonian Hs is not 
bounded below, and there is no ground superstate. The thermal situation does not allow 
the asymptotic particles. The requirement of the asymptotic particles at finite tempera- 
ture leads to the trivial S-matrix 111 0]. Many works on thermal field theories respect the 
Heisenberg equation for the field operators, but do not make the representation space 
clear. We will provide a definite representation space in the Liouville space, constructing 
the unperturbed representation there, which follows from the very basic requirements. 

Let us put the following three requirements. 

(a) The unperturbed representation of the superoperator formalism is constructed on 

the quasiparticle picture at each instant of time. 

(b) In the causal Green's function, the macroscopic time-dependent quantities affect 

the microscopic motions only in the future but not in the past, which we call the 
thermal causality. 

(c) The system approaches to equilibrium after a long time, consistently with the law 

of thermodynamics. 

The requirement (a) realizes the concept of the quasiparticle by Landau in our treat- 
ment. We restrict ourselves in this paper to the case in which the original Hamiltonian 
possess the global phase symmetry, and its spontaneous breakdown or its rearrangement 
does not occur. (The present construction of the unperturbed representation will be ex- 
tended to cases with spontaneous breakdown of symmetry or symmetry rearrangement in 
the future work. ) Consider the unperturbed system in the Schrodinger picture. Then (a) 
implies that the unperturbed superstate |po(^)} exists for any t . The superstate |po(i)^ 
which we seek is a special one, playing a role of "vacuum" in the Liouville space on which 
the quasiparticle picture is constructed. As the global phase symmetry of the system is 
kept, we require that |/f>o(0)) should be invariant under the global phase transformation, 

cij e l9 a,j , a,j -> e~ l9 a , (44) 
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and their Hermitian conjugates. Then its general form is 

m, m) . (45) 

m 

The three basic properties of the density matrix are (i) the Hermiticity pj(i) = po(t), (h) 
the normalization Tr[po(^)] = 1 and (iii) the positivity (</>|/?o(^)| ( / > ) > for any \<fi) , and 
their expressions in terms of p m (t) are 

p* m {t)=p m {t) for(i), (46) 
Y,Vm{t) = l for(ii), (47) 

m 

Pm(t) > for (iii) . (48) 

The mean unperturbed number nj (t) is given by 

nj(t) = Tr[a]a jPo (t)} = (l\a]aj\p (t)} = £ m jPm (t) . (49) 

m 

The unperturbed super-Hamiltonian, which is diagonal for each mode, at most bilin- 
ear in cij , alj , dj and , and invariant under the transformation in Eq. f)44[) , is generally 
expressed by 

= ^ ^rjxj{t)a^/aajdj + r)2j{t)a^/aa\d\ + //:;,(/ + ^-(^ajaj + V5j{t)} , (50) 

3 

where the operators are time-independent (the suffix S is suppressed) but H u (a,--- : 
r)(t)) becomes time-dependent through the five undetermined parameters riij(t) (i — 
1,2,- •• ,5) for each mode, complex in general. Note that H u (a,--- : i](t)) is non- 
Hermitian. The requirement (a) also requires the unperturbed Liouville-von Neumann 
equation, 

i^\p (t)}=H u (a,---:v(t))\po(t)}- (51) 
For the property (i) to be consistent with Eq. (f5"Tj) . one has 

{H u (a, T)(t))}~ = -H u (a, f](t)) , (52) 



as Eq. ((22]) and the tilde conjugation of Eq. (fSTj) show. The condition (|52l) implies that 
the five complex 77y(i)'s are expressed by five real parameters for each mode 

{Wj(*),Cw(*),0«(*).C3i(t),CBi(«)}, 

Vij (t) = i&j (t) , mj (t) = i&j (t) , 

V3j (t) = Uj (t) + i( 3 j (*) , Tj4j (t) = -UJj (t) + i( 3j (t) , 

V5 j (t)=<5j(t). (53) 
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The conservation of probability (the property (ii)) is guaranteed if 

(l\H u (a,---:r)(t))=0, 
which is seen from (|5ip . From this follow the two constraints 

C« (*) + <T&i (*) + 2(3j (*) = , & (t) + C 5j (t) = . 

Using Eqs. (|49p and (j5"Tj) . one can calculate nj(i) as 



(54) 



(55) 



^ -(l\a]a 3 H u (a, • • • : »/(t))|po(t)> - -Cy(*)»,-(i) + &(*) (1 + «*,•(*)) • (56) 

Equations f|55[) and (1561) form a simultaneous equations for the four unknowns d(t) 
(i = 1, 2, 3, 5) and their general solutions are 



CyW = <™j(*)+7i(*). 

1 + 2cTOj (t) 



•;(*) 



^(*) = -0y(*) = ^) + I+t^w 



C3j(*) = ~ohj{t) 



2(1 + <rnj(t)) 



7i(*) ; 



(t)7i(*). 



(57) 



with an arbitrary function 7j(t) . The parameter is also arbitrary at this stage. 

Let us move to the Heisenberg picture for the unperturbed system, or to the inter- 



action picture with the time- independent superstate |po) 
equation, equivalent to Eq. (|51l) . is 



|po(to)} ■ The Heisenberg 



i-A[t) = [A(t),H n (t) 



A{t) = U^{t,t )A s U {tM) 



(58) 



where the definition (1321) has been used and 



i 

+ i [Cy {t)<ry/ahi (t)5j (t) + Caj (t)<^aj (*)ot (t) 

+Cm(*) (aj(t)aj(*) + aj(t)a i (t)) +Gy(*)} 



(59) 



with Eq. ([5"7ji . 

Although the availability of the Feynman diagram method in calculating the causal 
Green's functions seems a technical matter, it is crucial. Without the Feynman diagram 
method, a systematic analysis of the Green's function would be impossible, since we 
could not use the powerful tools of quantum field theory such as the Wick theorem, 
the Dyson-Schwinger equation, the concepts of one-particle irreducible diagram and the 
effective action, the renormalization procedure and so on. For the availability of the 
Feynman diagram, one should have oc (l\ in Eq. (|32]t. This relation holds if the 

interaction super-Hamiltonian in Eq. (|39p satisfies 



[l\H 1 (t) = 0. 
10 



(60) 



As it can be shown that ^I\H(t) = for the Hermitian H, we only have to check 
(l\H u (t) = 0. Equation ([54]), coming from the conservation of probability, is true for 
arbitrary nj(t) , 7j(i) and u)j(t), and H u (t) is related to H u (as, ■ ■ ■ : r)(t)) as in Eq. (j35|) 
with Eq. (|32|) . so we have 

(l\H u (t)=0. (61) 

Thus the Feynman diagram method is available. 

The mean number rij(t) is a macroscopic quantity, given by the expectation over 
the density matrix po(t) . It is expected that our unperturbed Hamiltonian governs the 
microscopic quantum dynamics in a consistent manner with the macroscopic thermal 
change in the background. The dependence of H n (t) on rij(t) means that the operators 
in the interaction picture {a,j(t) , dj(t) , 2j(i) , oJj(t)} depend on rij(t) . Suppose the two- 
point unperturbed causal Green's function or the unperturbed propagator, 

A(*i,ta) = -i(/|T[ai(ti)a 2 (*2)] \po) 

= -id(h - *2)(/|oi(t 1 )a 2 (t2)|po> ~ ivBifr - ii)(/|a 2 (t 2 )ai(t 1 )|/Jo» , (62) 

where a, stands for any of {hj{t) , aj(t) , 5j(t) , Oj(t)} . Generally speaking, A(ii,t 2 ) 
depends on t\ and i 2 , and on both of ny(ii) and %(i 2 ) ■ The requirement (b) demands 
that it should depend only on nj(i 2 ) but not on nj(t\) for ti > t 2 and vice verse for 
t 2 > t\ . Otherwise the microscopic dynamics would be affected by the macroscopic 
quantities in the future. 

Let us see what constraint emerges when (b) is imposed. We have the three parame- 
ters uij (t) , rij (t) , 7j (t) . The parameter ujj (t) is interpreted as the renormalized excita- 
tion energy and belongs to the microscopic dynamics. On the other hand, rij(t) and "fj(t), 
appearing in the imaginary part of the super-Hamiltonian, represent thermal situation, 
and therefore are macroscopic parameters though the physical meaning of 7j (t) is not 
clear at this point. To extract the macroscopic part in the unperturbed Hamiltonian, we 
introduce new operators {<5y(t) , &j(t) , Oj(t) , a^(t)} from {&j(t) , &j(t) , a,j(t) , , 



The unperturbed super-Hamiltonian for the new operators is 



(63) 



H au (t) =iJ2 [Ci 3 {t)a^a 3 {t)& 3 {t) + C, 2j {t)ay/Za](t)a)(t) 

3 

(tf{t)&{t) + a\{t)&i{t)) +CsiW 
with Eq. (I57[) . and their equal-time (anti)-commutation relations are 



(i) , a£ (t) = 5jfc , others = . 



(64) 



(65) 



Replace a,i(t) with a!i(t) in Eq. (H)2"|) and consider for ti > f 2 . As the time-dependence 
of &i(t) comes solely from rij(t) and 7j(i) , it is sufficient for (b) that ^/|&i(ti) is inde- 
pendent of t\ . Similarly for £ 2 > t\ ^/|ci! 2 (t 2 ) is independent of £ 2 . Thus the requirement 
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(b) is fulfilled if we have all of the following relations: 

= = (I\&](t) = (I\&j(t) = (I\a](t) . (66) 

The equations for a(t)'s are 

&j(t) = Oy(t)trv^fi}(*) + di(t)&j(t) , 

= -&,■(«) v^5 3 -(*) - Cy (*)&}(*) , 
= C?i(t)v^<sJ(*) + C« (*)*,■(*) , 
°j (*) = -Cij (t)ffv^<*j (i) - C 3j (t)fij (i) , (67) 
from the Heisenberg equations i&j{t) — \&j{t) , H an (t)] and so on. We also note 

(/| (&j(t) - a^&){tj) = (l\ (aj(t) - aV^fij(*)) = , (68) 

derived from Eq. (JTJjJ and (l|iT au (*) = . We find from Eqs. (J6Z) and (65) that Eq. ([55]) 
holds true only if cr^j (i) + (3j (t) = and £ij (t) + C3j (t) = or equivalently 

7,(0=0. (69) 

This way the unperturbed representation can be described by only the two physical 
parameters ujj(t) and rij(t), which reflect thermal situations of the system and can be ad- 
justed freely to them. Explicitly the unperturbed Hamiltonian in the interaction picture 
is 

H u (t) = [wj(t) - 

3 

+ihj(t) U/aaj^ajit) + cxy/a&X{t)a]{t) - oa]{ty aj {t) - a 3 -(f)at(t) j , (70) 

or 

H au (t) = i ^2 n>j(t) {v^&j + cr^/aaX(t)aX(t) — aaX(t)dij(t) — Oj(t)aj(i) j 



5. Derivation of geometric distribution 

We study the unperturbed dynamics obtained in the previous section. 
Substituting Eq. (|69l) into Eq. ([57]) and rewriting Eq. (|67p . we have 



&j{t) — —ahj(t) (aj(t) - a\/aa]j(t)^ , 

aj(i) = -a^hj(t) {&j(t) ~ a^a]{t)\ . (72) 
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It is easy to see that 
d 



(73) 



This can be understood from the fact that dtj{t) — *Jaa]j(t) and its tilde conjugate com 
mute with H au (t) in Eq. (|7T|) . 

We also have the following time-independent combinations, 



d_ 

dt 



{(1 + oMt))&i(f) - a^n 3 {t)a]{t)} = - {(1 + ^ (t) - \/^ (*)<#)} = . 

(74) 

The proof of the first one is as follows: 



- {(1 + on, (*))*,(*) - (*)«}(*)} = <m*(t)- {fy(t) - 

+ (i) + crnj (t) |aj (i) - v^5] = , (75) 



due to Eqs. (f?2"j) and (j?3")) . The second one can be proven similarly. 

We shall elicit the implications of Eq. ([74]) for boson (c = 1). For simplicity, we here 
discuss a model of a single mode, but the result can readily be generalized to multi-mode 
cases. Noting that 



= U a ^(t,t )aU a o(t,t ) , 



U a0 (t,t )=T 



exp — i / H au (a, • • • : n(s)) ds 



and 



(76) 
(77) 

(78) 



\p (t)} = J7 a o(Mo)|po) = y^p m ft)[m, m) . 

Operate {(l + n(t))a(i) — n(t)er (i)} on \po}, then we have 

{(l + n(t))a(t)-n{t)&\t)} |p » 

= fr-^i, t ) ^ vWl { (1 + n{t)) Pm+1 {t) - n(t) Pm (t)} \m + 1, m» . (79) 



Equation (|74j) implies for arbitrary ii and 

^2q m (ti)\m+ l,m) = y^gm(*2)#a(fri,*2)lra + 1, 



g m (t) = Vm+l{(l + n(t))p m+ i(i) - n(t)p m (*)} (m = 0, 1, 2, • • • ) . 
Putting t\ = t + At , £2 = t and taking the limit At —> 0, we derive 



(80) 
(81) 



q m (t) = h(t)Vm + l{Vm + 2 <7m+i(i) - 2\/m + 1 g m (*) + \/m g m _i(*)} 

(m = 0,1,2, •••), 
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(82) 



or equivalently 

q(t) = n(t)Mq(t) , (83) 

with a real symmetric matrix M 

{-2(m+l) (m = to') 

^(to + 1) (to' + 1) (to = to'±1) (to = 0,1,2,---)' (84) 
(otherwise) 

Any q at arbitrary time t can be expanded as 

q(t)=J2 c t{t)ui, (85) 

with the orthonormal-complete set of eigenvectors, {ui} : 

Mui = \tUg . (86) 

Solving Eq. (j83|) . we obtain 

^^qM^^N-w] , (87) 

i 

which implies 

\\q(t)\\ 2 = Y,\ C ^ T )\ 2 e2X ^ n{T) ^ ■ ( 88 ) 

£ 

The trivial solution of Eq. (JH3J) is 

q(t)=0, (89) 

or equivalently = for all £, which in turn leads to the geometrical distribution for 

Pm(t) 

= TT^TTT = f® (^dependent of to) , 
Pm(0 1 + n(t) 

Pm(t) = (l-f(t))f m (t). (90) 

Equation (|88|) indicates that if q(r) ^ at a certain time r then q(oo) 7^ 0. If we 
require the approach to equilibrium (the requirement (c)) for which q eq = , then only 
the meaningful solution is Eq. (|89[) and p m (t) is the geometrical distribution, 

\ Po (t)} = (l-f(t))Y / f m (t)\^™}- (91) 

rn 

For the fermionic case, there are only po(t) and so we always write 

|po(t)» = (l + /(t)r 1 E/ m (*)l m ' TO »- ( 92 ) 

tn 
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Equations (|9T|) and (|92|) are unified into a single expression, 

lM*)» = (i-^/(«)rE/ m wi m ' m »- 



(93) 



Thus for the time- independent superstate |po)), consistent with the requirement (c) 
we have 



{(1 + anj{t))6ij{t) - o-^nj(t)at(«)} \po) 
{(1 + o nj {t))a 3 (t) - V^i(t)&t(t)} |po» = . 



(94) 



6. Equivalence to nonequilibrium Thermo Field Dynamics 

We rewrite the superoperator formalism, obtained from the basic requirements, in 
language of Thermo Field Dynamics The operators {dj(t) , Oj(t)} and {dj(t) , Oj(t)} 
in the interaction picture can be identified as the non-tilde and tilde ones in the nonequi- 
librium TFD formalism, which form the thermal doublets, 



aS(t) 



d 3 (t) 



a*(t)= d](t) -y/Ba^t) , (M=l,2) 



with the equal-time (anti)-commutation relations, 

[a£(i) , a%(t)] a = S jk S MV , [a£(t) , jJ(i)] ff = [a£(t) , aftt)], = . 

The superstates and \po~} become the thermal vacua, denoted simply by (i 
for which the following thermal state conditions hold, 

(0| (aj(t) - a^a){t)) = (0| (a){t) - a^a^t)) = , 

{(1 + anjitycijit) - (T^^^t)} |0) 

= {(1 + anjWajit) - ^n 3 {t)a){t)) |0) = . 



(95) 

(96) 
and |0), 

(97) 
(98) 



from Eqs. JBSJ and 

The properties in Eqs. ((97)) and (l98l) enable us to introduce the time-independent 
annihilation- and creation-operators associated with the thermal vacua, , , ^ , <fj} 

&|0>=&|0>=0, (0||] = (0||]. 



(99) 



with 



(*) = 



.V^}(i). 



= Bf(t)o^(t), ff(t) = ^(t) -v^OW =^(t)B7 1 ^w, 



(100) 
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where 



&(<) 



-i J t uJj (s) ds 



£,•(*) 



-i f t Wj (s) ds 



— Jj tuj j(s) ds 



(101) 



Note that the dummy thermal index is summed in what follows. The thermal Bogoliubov 
matrix B^ v (t) is defined by 



Bf (t) = 



1 + arij{t) —crrijit) 
-1 1 



and the (anti)-commutation relations of the ^-operators are 

. = *i* v . k; , ^ = , si* = o . 

Equation f|100[) can be inverted into 

= Bj^{t)%{t) , ojf(*) = ejwsf (t) , 

which are the starting relations in nonequilibrium TFD Q. 

The unperturbed and interaction Hamiltonians are Eq. (|70[) and 

#i(t) = # - # u (f) 

H = H + H[ nt 

H a = J2"aj {a](*)4i(t) - a}(*)fii(t)} , 



(102) 

(103) 
(104) 



(105) 



respectively, where 75/i nt = i?i n t ~ ffint represents a non-linear interaction. Note that the 
counter terms, proportional to 8uij(t) = u)j(t) — Uoj and hj(t) , are present in Hi and 
that they will be determined from the self-consistent renormalization condition. In the 
doublet notation, H u (t) is 



(106) 



where 



Tn = 



The causal Green's function for the Heisenberg operators is defined by 



G£(ii,t 2 ) 



< H (*iKh(* 2 ) |0), 



and can be rewritten in the interaction picture as 



Sa^h)aUt 2 ) |0) 



(107) 



(108) 



(109) 



because of <0|5'^ 1 = (0| due to Eq. {60]), where S was defined in Eq. (03) with Eq. (l40l) . 
Because of the Dyson's expansion formula in Eq. (| 109[) and the existence of the £- 
operators, the Wick theorem holds. Therefore the Feynman diagram method is available 
in calculating the causal Green's functions in nonequilibrium TFD. 
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7. Self-consistent renormalization condition 

According to the Feynman diagram method, we have the Dyson-Schwinger equation 
for the two-point Green's function, 

/OO />OC 
OO J — OO 

J ) K 

(110) 

where the unperturbed propagator is 

A£(ti,t 2 ) = -<(0|T K(^K(t 2 )] |0) . (Ill) 

We also define the unperturbed and full propagators for the ^-operators by 

d%{tiM) = -*<0|T[#(t 1 )g(t 2 )]|0) , (112) 
9%{t x ,t 2 ) = -i(0\T[$ H (h)% H (t2)]\0) , (113) 

respectively. These propagators are related to those for the a-operators as 

A£(ti,t 2 ) - Br 1 >w'( tl )dfc'(t 1 ,t 2 )B u k ' u (t2) , (H4) 
G^(t 1 ,t 2 )=BT^\ tl )g^'(t l ,t 2 )B^(t 2 ). (115) 

While the unperturbed propagator d has a diagonal structure with respect to the thermal 
index, 

e -*f£ (116) 

the full propagator g has an upper triangular structure in general, that is, gjl(ti,t 2 ) cx 
9(ti-t 2 ),g^(h-t 2 ) cx e{t 2 -t 1 ),g]l{t 1 ,t 2 ) ^ and gf k {t u t 2 ) = 4]. This is because 
£j H and £j H identically annihilate the bra- vacuum, see Eqs. (l60l) and (|68|) . It follows 
from the tilde-conjugation that has the following properties: 



9-l(h,t 2 ) = glf*(t 2 ,h) , gf k {t x ,t 2 ) = -9%'*(t2,h) ■ (H7) 

The self-energies, defined through the Dyson-Schwinger equations S = A -1 — G _1 
and S = d^ 1 — g^ 1 , have the following properties, 

^(t 1 M=BJ 1 ^\t)S^'(t 1 ,t 2 )B^(t'), (118) 

[ S^(ti,t2)J ' ( } 

S#(fi,f 2 )oc0(ti-f 2 ), S$(ii,t 2 )oc0(t 2 -ii), (120) 

^(ti,t 2 ) = Sg'*(t 2 ,ti), Sj|(ii,t 2 ) - -Sg'*(t 2l i!). (121) 



We move to perform the self-consistent renormalization, that is, to fix the two pa- 
rameters u>j(t) and rij(t) in the unperturbed Hamiltonian H u (t). 
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d%(h,t 2 ) = 6 jk 



-i9(h - t a ) 
iOih-ti) 



First recall the on-shell renormalization condition in equilibrium case which is a stan- 
dard prescription in ordinary quantum field theory and fixes the renormalized energy ojj . 
Then the self-energy depending only on t\ —t2, S^(ti —1%) , is given in fco-representation 
by 

/oo 
drS%(T)e> k °\ (122) 
-oo 

and we have the spectral representation in equilibrium case which becomes in t- 
representation, 



/oo 
dn B- 1 ["(«)]""' 
-oo 



^0 (T) i9(-r)f V BM*)f v (123) 



with the equilibrium distribution n(n) = \j(eP K — er) for the inverse temperature /3, 
and the spectral function ajk(n). The j-diagonal elements in fco-representation have the 
properties of 

S£(ko) = S%'*(ko), (124) 
Sjf(k Q ) = 2iaTr(n(k ) - n 3 )a 3] {k ) , (125) 

and identically Sjj(k ) — as in Eq. (|119D . Because nj is the equilibrium distribution, 
Sjp(ko) automatically vanishes at on-shell: fco = ujj. Finally we require the on-shell 
renormalization condition, 

/oo 
dr Sjj(T)e l ^ T , (126) 
-oo 

to determine uij. Equation f|126[) implies the condition Re Sj*(ujj) — as well because 
of Eq. (|124p . This way all the elements of the self-energy Sjj (uij) become zero, except 
for the imaginary parts of Sjj(ujj) and Sjj(ojj) representing thermal instabilities. 

The self-consistent renormalization condition in the equilibrium case can not be ex- 
tended straightforwardly to the nonequilibrium case. Firstly, since the time-dependence 
of the self-energy is not simply t\ — ti but it depends on t\ and t% separately, its ko- 
representation and the definition of the on-shell are not trivial. Secondly, an additional 
condition is required to determine hj (t) . For this Chu and Umezawa has proposed their 
renormalization condition not on the self-energy but on the propagator, called diagonal- 
ization condition [Bl. Il9l. 18 1. The diagonalization condition is 



$#(*,*)= 0, (127) 

and determines the temporal evolution of rij(t) in the leading order. To show it, we solve 
formally the Dyson-Schwinger equation for g 12 , and obtain 

9jk(h,h) =Y1 [dsid 2 g}},(t 1 , Sl )S}? kl ( Sl , S2 )gf? k (s 2 ,t 2 ) . (128) 

Next, we approximate the full propagators g 11 ' 22 - ) in Eq. (|128[) by the unperturbed ones 
d 11 ' 22 ) and divide the self-energy S into a loop contribution Sioop and a contribution of 
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the thermal counter term Sq, 

^ifc,Q(*i>*2) = -io-hj(ti)5jk$(ti - t 2 ) 

Then we have 

gfAt,t) = -i 



1 




i if 



ds 



ahj(s) - 21m / ds S jjAoop (s, s ) e 



'1 P i Ss' ds " w i( s ") 



(129) 



(130) 



and the diagonalization condition Eq. f|127[) implies the following transport equation 



fij (t) = 2aTm / ds S 



? 12 



(i,s)e 



i J t ds' ujj (s') 



(131) 



which determines the temporal evolution of rij(t). The loop contribution Sjj loop can be 
calculated diagrammatically. 

Although the transport equation, obtained by the diagonalization condition (|127l) . is 
reduced in the Markovian limit to the ordinary quantum Boltzmann equation derived in 
the other methods, there is a critical problem, that is, the diagonalization condition in 
higher orders leads to be inconsistent with the equilibrium theory. To confirm it, we ex- 
press gjj(t, t) in terms of the unperturbed number distribution rij(t) and the Heisenberg 
number distribution nju(t). The latter is defined by 

n jH (t) = (Q\a} H (t)a jH (t)\0) , (132) 

or equivalcntly 

G)){t,t + 0) = -ian jH {t), G)){t,t - 0) = -t(l + an jH {t)) . (133) 



On the other hand, we obtain from Eq. (|115[) 

G)) (t, t + 0) = -ian 3 (t) - g$ (t, t + 0), G)} (t, t - 0) = + artj (t)) - c/j| (t, t-0), 
and therefore 



(134) 



■(t,t)=ia(n jH (t)- nj (t)) 



(135) 



Thus, if the diagonalization condition is applied, the Heisenberg number distribution has 
to be equal to the unperturbed one. But it is not true in the equilibrium limit, where 
the full propagator has the spectral representation, 



with 



dn pjk{n) = 8jk and 



dn 



B-'Hn)} P JM.. B[n( K )] 

Kq — K + IET3 



(136) 



dn n(n)pjj{K) , 

o 
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(137) 



and we have rijH ^ n(uij) in general. 

We propose a new self-consistent renormalization condition here. For this, the on- 
shell energy renormalization condition of the equilibrium Eq. (|125|) has to be extended 
to the nonequilibrium case, and the diagonalization condition Eq. (|127l) has to be re- 
placed with a new one which yields the correct equilibrium limit. It is desirable from 
a unified viewpoint that the two conditions have similar expressions, in contrast to the 
situation that the diagonalization condition is imposed on the propagator while the en- 
ergy renormalization is the condition on the self-energy. We attempt conditions on the 
self-energy. 

A natural extension of the on-shell self-energy as in Eq. (|122p with ko = ojj to the 
nonequilibrium case with the time-dependent uij (t) would be 

/oo t 
drS^{t 1 M)e i ^ dsUi{s) , 
-oo 

/oo t 
dr {e(T)S^ r (h,t 2 ) + e(-T)S^ a (tx,t 2 )} e'/i^W (138) 
-oo 

with t = ti — t 2 . There is an ambiguity in this expression, namely, as to how t\ and t 2 
depend on t and r. Let us parameterize 

t l= t + r t T, t 2 =t-(l-r e )r (l = r,a), (139) 

for retarded and advanced parts separately, where ri are undetermined parameters. We 
expect that 0Jj(t) and fij(t) are determined from S^[ojj\£\, and the integrands in (|138[) 
should not involve the quantities later than t from the viewpoint of the thermal causality 
as was already discussed, which implies that t > t\ , t 2 and therefore > rgr and > 
— (1 — ri)T. We have the consistent choice, r r — for the retarded case r > and r a = 1 
for the advanced case r < 0. Thus we acquire the unique expression 

/oo 
dr ld(j)S% r Qt,t-T)e i f*-T d ' u 'W +9(-T)S%Jt + T,tyti +Tds ^ s n . (140) 
-oo 

Now we propose here a new set of renormalization condition, a successor to the diago- 
nalization condition (|127[) . as 

Sj? [lu 3 ;t]=0, and Re S]] [lu 3 ;t]=0. (141) 

With the aid of Eq. (|12ip . we obtain 

ft 

S}f[LOf,t] = 2Tm / ds S^(t, s yf° ds ' ^ (s,) , (142) 

J — oo 

which with the on-shell renormalization condition (|141[) implies the same transport equa- 
tion as Eq. (|13ip . Note however that it never means the two conditions (I127p and (|141j) 
are equivalent. The leading order replacements g 11 ^ 22 ) — > d 11 ^ 22 ^ in the left-hand side 
of Eq. (|128p were necessary to derive the transport equation under the diagonalization 
condition. While the diagonalization condition gives a correct transport equation only 
in the leading order, our new on-shell renormalization condition is valid in any order. 
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8. Summary and discussions 

In this paper it was shown in the superoperator formalism that the three basic re- 
quirements on thermal field theory restrict the structure of unperturbed representation 
in the interaction picture uniquely, that is, the corresponding unperturbed density ma- 
trix has the geometrical number distribution Eq. (|93[) . From this representation follows 
nonequilibrium TFD using the thermal vacuum and time-dependent Bogoliubov matrix, 
though it has been an assumption in the nonequilibrium TFD. 

The three basic requirements are (a) the existence of quasiparticle picture at each in- 
stant of time, (b) the thermal causality that the macroscopic quantities should affect the 
microscopic motions only in the future but not in the past, (c) the relaxation to the equi- 
librium after a long time. While the requirement (c) is the thermodynamical law, (b) is 
necessary to provide a consistent causal description in the coexistence of microscopic and 
macroscopic quantities. The representation space or the particle picture for the nonequi- 
librium system of quantum field with change in the background number distribution is 
made clear due to the requirement (a). From the viewpoint of the loop calculation, the 
use of the time-dependent unperturbed superstate (or the unperturbed density matrix) 
provides the better approximation than that of time-independent one, usually fixed at 
the initial time, because the former takes account of the temporal change in the number 
distribution. In connection with (a), we point out that the unperturbed Liouville-von 
Neumann equation (j5Tj) with Eq. (|70|) is a special form of the Lindblad equation 20]. 
The Lindblad equation is derived for a system interacting with an environment system 
after some coarse-graining in the time direction is performed. We consider that following 
our time-dependent unperturbed density matrix involves coarse-graining process. 

Since the unperturbed representation is written in terms of the two unspecified pa- 
rameters, the unperturbed number distribution rij(t) and the excitation energy u>j{t), 
the expectations of the Heisenberg operators can not be obtained at this stage. The way 
to fix the parameters is known and simple in the equilibrium case, namely, rij is nothing 
but the equilibrium distribution l/(e^ — a), and ujj is interpreted as the renormal- 
ized excitation energy and is determined through the on-shell renormalization condition 
(|126p . It is non-trivial to extend the definition of the on-shell in the nonequilibrium 
case where the time translational invariance is broken. Chu and Umezawa proposed the 
diagonalization condition which was not on the self-energy but to on the propagator, and 
derived the quantum transport equation in the leading order. They have not given an 
explicit method to determine uij(t), and it was approximated by the bare energy. We 
have pointed out in this paper that the diagonalization condition becomes inconsistent 
in the equilibrium limit in higher orders and have proposed a new renormalization con- 
dition which is a natural extension of the on-shell renormalization in the vacuum and 
equilibrium theory and determine rij(t) as well as ujj(t) in any higher order. The concept 
of the thermal causality removes the ambiguity in defining the on-shell in nonequilibrium 
situation. 

Our discussions throughout this paper are restricted to the cases where there is no 
spontaneous breakdown of symmetry. When a symmetry is spontaneously broken and 
a Bose-Einstein condensate is formed, we introduce the Bogoliubov transformation to 
diagonalize the unperturbed Hamiltonian having terms such dj&j. There are various in- 
triguing phenomena in the cold atomic gas systems associated with the broken symmetry 
typified by the unstable condensate decay and the quantum phase transition. We know 
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the two methods to deal with the spontaneous breakdown of symmetry in nonequilib- 
rium TFD. The first method is to introduce a 4 x 4- matrix transformation to combine 



the thermal and usual 2 x 2-transformations [21| . But then there are so many parameters 



in the 4 x 4-matrix that we have not succeeded in extracting physical parameters yet, 
as we have extracted the two physical parameters Uj(t) and rij(t) in the 2 x 2-matrix. 
The second method is to expand the field operator in the complete set of time-dependent 
wave functions affected by the time-dependent order parameter [2"3.Hoj. In this method, 
the 4 x 4-matrix is given as a direct product of the two 2 x 2-matrices. Both of the 
methods are not yet derived from the superoperator formalism, and the extension of the 
discussions in this paper to the symmetry broken cases is a future task. 
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